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ABSTRACT 


If a weapon system is fired at a circular target and the 
impact point is distributed as a bivariate normal random 
Marta le, 26 15 not particularly difficult to determine hit 
probability if the expected impact point is at target center. 
If, for some reason, the expected impact point is offset, the 
problem of determining hit probabilities becomes quite complex. 
Herein is developed a method to approximate such a hit 


probability . 
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eee RODUCT LON 


Pom oeObaomt ty Of wmaitting a Specified target with a 
given weapon system is of particular interest to both mili- 
baby Commanders and weapon systems analysts. Determining hit 
probability is not always a Simple matter. In particular, if 
the impact points of the weapon system are distributed accord- 
ing to some probability distribution, and the expected impact 
Boimt iS not at target center, the pEooien of determining hit 
probabilities becomes quite complex. 

Over the past thirty years this problem has been studied 
in depth as the extensive bibliographies of F. Grubbs [Ref. 
ieee enek ler [Kef. 2], and W. Guenther and P. Terragno [Ref. 
3] illustrate. Several approaches to solving the problem have 
been investigated by these and other authors in an effort to 
Simplify the computations required to estimate or actually 
compute the hit probabilities. 

As an example of one approach to the problem, Grubbs [Ref. 
1) uses what he calls a “Single, straight-forward and rather 
Simple technique" to approximate these hit probabilities in 
which he employs an approximate central chi-square distribu- 
tion with fractional degree of freedom, or a transformation to 
approximate normality. There are, of course, other methods 
available to approximate these hit probabilities, one of which 


will be discussed in the following pages. 
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The method which will be discussed is computationally 
easy, using simple mathematics in conjunction with a table of 
modified Bessel functions. These tables are usually avail- 
able to the analyst and enable him to accurately approximate 
the desired hit probabilities without the use of complex com- 
putations. It should be mentioned that tables of offset 
circle probabilities have been compiled by Rand Corporation 
[Ref. 4]. These tables give hit probabilities as a function 
of target radius and offset distance of expected impact point, 
but are not as readily available Berchic Bessel function tables 
used in the following procedure. It should also be noted 
that routines for Bessel functions are available for most all 


Pomeoucet Systems. 





II. NATURE OF THE PROBLEM 


Semele tecamecrrecilak target with radius a. Place the 


Paioin OF ene coordinate system at the center of the target. 





Let f(x,y) be the probability density function of the 
gmpact point, (x,y), of the projectile fired at the target. 
Let (my, My) be the mean of the impact point distribution and 
oR and Ty be Eiomvartanece or the 2mpact point distribution 
in the x and y directions respectively. Assume covariance is 
Zero. tf the distribution Of impact points follows a bivari- 
decemormal distribution, then 


| ans v-tay\~ “- : 
Fly = = ene (-a[( a (8) 
LOT, 2 
a Zz 2 
PosuNew2— = OF , and R= My + My ees oe be the proba- 


Daeheyeewedeitit On thaevtarget, then 


= Ze 
a ork } if exe boa |( eo) (y-Hy) } dy ely 
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Changing to polar coordinates yields 
t= Gocnee y= r sin® 


and 


Zh 6 
he — Jeo) exe} J a[(recs@ mys (r-6in 8 -uy\]} dr 
ER: Lr 
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Rotate axes so that m,>= R and By = O , then 


ge Aa z 2h 


Pp. eo" ( eve(%)) feve ("BE Bese) aol dr 
oO 
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Heetong at the term, fexe(tR jae 
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Lh 
and {cos 0d = (="] cos Bab 
a) ran 9 
ein 1S odd, tnen [eos” GAB =O . If n is even, then 
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pe 94 0-iactions er ie 
(MN -2)(N-4) 4-2 
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4K ~ hus; 


> 2h a ry" Fenty? 2 Also, for n even and equal to 2m, 
(K-15) FT 2g es LNA = Meee 
nl(n-2\n-4)-°°4-2 Tyv(n-2yin-4): -4.2)° (yes 
(2m)! 2% 
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The Bessel function, denoted by Jp (¥) Meas Legual co 


ad nh Ys ad 
3 (“0 iz) , where x is the argument of the 
> ni Msp)! 


Bessel function and p is the order of the Bessel function. 


Miewmodiried Bessel function of order 2eNeC is denoted by ie (x), 


and is equal to un (tx) = 57-0" (8 () - <{ Gy 


Nzo —e N=O T™miy* 


2 ir 
thus, feep (cos 8)cld = nt (1 “yt = 2% T,¢¥), 


Ta 
so that feve("Gaces O\d@ =Z2% T, (8) 


v 
It was shown previously that the peop ley et Hitting. ene 


ee 
target, P , Was equal to Gay r* Srevel aba [crete ‘Ja dlr. 
ZnO 
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Lh 
iceewas also shown that feve lS Bay oo « rod ial Be (18 ). 
an 


Bs a result, P = em [rete Ze (% 8 \dr or 
znor J 
als yr 
p= e2o" |® en ie (8 \dw 
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It 1s only necessary to know the ratios, = and = ao 


determine the probability of “ok EnewEarge:, 
-~ pt ye 
Rw a gt 27+ 
Plgiel? <s fre ae VS S)dr. 
o 
The Rand Corporation has constructed tables giving these 
hit probabilities as a function of the two ratios, © and = ; 
T 


ier. 3) . 


Gilliland, [Ref. 4], shows that 


Se es tas (&,.) 


where ie Cae =j- ro (a i} a 


N=0 
make, ([Ref. 5], Pee the function 
\ y 
J(u,4) = 1- 33 a 
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~-U 
.2¢ = (Sag 
EST, JH Ee)= {- 2. e fre a c* Ie (S42) dae 
As a result, 


eo <a 
Ro) = 41-J(#. 8). 


Luke, [Ref. 5], gives the expansion of J (4,4) ina 


eemtes of Bessel functions. 





Sly) = en T, (8) for 24 
J (4,4) = 1- as oat r,, (6) for Y >1 
where G = Zly4)> ; and " e (ay 
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These latter two equations yield the true hit probability 
elec Weapon System fired at a circular target with radius a, 
and normally distributed impact points with expected impact 
point at Py and Hy ° 

A good approximation to the true hit probability can be 
had by considering only the first few terms of the infinite 


series in the appropriate equation. That is 


Cee 


a < Q 4 
8 (8,$\=1-€ = Nal Ly (752) Rs 
: (38) oA avr ra R R >on 
7 € 2) aS | 


R21 


where VY\ iS, hopefully, relatively small. If WN is less 
than five, then the approximation would be beneficial to 
analysts and would be quite simple to use. If WM is 

greater than five, the calculations required to get a reason- 
ably accurate approximation to hit probability become 


SPronitbrtrve. 


dak 





fit eeONecLuslONS 


The accuracy of the above approximations was checked on 
an IBM 360 computer. The program used the FORTRAN language 
moomewoO LEM Subroutines, 10 and INUE, which computed the 
modified Bessel functions. This program is available through 
the Operations Analysis Department of the Naval Postgraduate 
Bemool in Monterey, California. 

The computer program was basically an iterative process 
in which different values of the ratios & ana 2 were set 
and the hit probabilities computed while the number of terms 
from the appropriate infinite series varied from one to ten. 
The approximate hit probabilities thus obtained were then 
compared to the true hit probabilities from the Rand tables 
[Ref. 4]. 

The two aspects of primary interest were (1) how accurate 
are the approximating equations in the various ranges of & 
and = , and (2) how many terms of the infinite series in the 
approximating equations must be considered before a reasonable 
degree of accuracy iS obtained. The specific values of the 
ratios = and = Mangement so tO 2,0 In the computer program. 
Botstrnenequation for R<a, the ratio a was initialized at the 
value .5 and the ratio = was allowed to vary between .5 and 
5.0. The value for = was then incremented to 1.0 and the 


process repeated with values of = greater than or equal to 


1-0. This procedure was continued until both e and a were 


a2 





Sgudeet© 5.0. For each value of . and : IS ShsyQicierayweuce 
hit probability was computed using only the first term of the 
infinite series initially, and then again for each additional 
term used, up to ten terms. In this manner it was possible 
to see how rapidly the approximating equations converged to 
eaes true hit probability. 

For the equation for R>a, the same type of procedure was 


used, holding the ratio 2 constant while the ratio ®& varied 


Ge rT 
through the appropriate values. As before, the ratio & was 
CT 


initialized at .5 and then incremented by .5 each time the 
ratio & Wad Gum its range of values and the hit probabilities 
computed. 

Representative values of hit probabilities obtained 
through the approximating equations are contained in Appendix 
A. Figure 1 compares approximate hit probabilities with ac- 
Etiabeehit PEObDabilities for R<a, and Fig. 2 does the same for 
R>a. In both cases the values shown were computed using the 
first five terms of the infinite series from the appropriate 
approximating equation. 

The convergence of the approximating equation values to 
the true hit probabilities is shown in the two examples be- 
low. Example one shows the convergence of the approximating 
equation for R<a as the number of terms from the infinite 


SSmiccmnereases. Values of & and O were selected to show 


_ 
both rapid convergence and relatively slow convergence. Ex- 
ample two shows convergence of the approximating equation for 


R>a. These examples along with Appendix A point out that the 





approximation TiMemoOBoodiotltty TOr R<a 1S Optimistic while 
the approximation of hit probability for R>a is conservative. 


Example 1 (R<a) 


eo So = 2.0 
o oe 
Approximation using: 

one term e200 69528 
two terms SELOS 707 
three terms L477 55 1 
four terms SLO 45.0001 
five terms 71044937 True probability: .104491 
six terms .- 1044937 
seven terms 10449 37 
eight terms , 044937 
nine terms Plo 44937 
ten terms 1044937 

Ro = 4.5,%= 5.0 

oe x 

Approximation using: 

Siew ce mm 9255550 
two terms -8596848 
three terms .8044760 
four terms ty CUES 
five terms » 7260436 lever wrooability: <653154 
Six terms TOs 0 Ss 
seven terms .6834279 
eight terms -O7 159.76 
nine terms .6639764 
ten Cerms Poa ga 7 le 
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Example 2 (R>a) 


Approximation using: 


one term .0690204 
two terms ROveazeo? 
three terms ~,0734669 
four terms 0734724 
five terms .0734725 PEleme rood lity: 2.073473 
Six terms 0734725 


seven terms GL BEI VA S| 
Sogitemeetuse 20734725 
Haine terms 40) 7) Be Ss 


ten terms "09754725 
oN 
= = 4,.5,8 = 5.0 
om om 
Approximation using: 
one term 0056702 
two terms L208 790 
three terms SG 525.04 
four terms PISS s713 
five terms 23304: Lane woroodo bility: | 2.272201 
six terms R419 2 70 


seven terms .2537574 
eight terms .2613786 
nine terms P2660 Go5 
ie (eerste = AS eee SIat 


Accuracy of the approximating equations can be improved 
in all cases by considering more terms of the appropriate in- 
finite series. If it is necessary to use more than five 
terms of the series, however, the mathematic manipulations 


involved become impractical to do by hand. Accuracy can also 


Ls 





oa a . 
be improved when or 1s equal to aot PEOMesOpendlx A tt 1S 


f 


apparent that when S 1s equal to = and both are greater 


eine s.,0, accuracy is the poorest. To more accurately ap- 
proximate the true hit. probability in this case it is suffi- 
cient to compute the arithmetic mean of the values obtained 
from both equations (R<a and R>a). For example, for 


R = oe — yy Ue 


\O 7 oO 
for R<a yields an approximation of .643, while the equation 


the true probability is .460. The equation 


for R>a yields an approximation of .325. The arithmetic mean 


of these two values is .484, which 1S more accurate than 


either of the two approximations. 
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IV. SUMMARY 


If the impact point of a weapon system is distributed as 
Peovariate normal random variable, the task of computing the 
hit probability for a circular target when the expected 


impact point is somewhere other than at target center is for- 





midable. The two approximating equations that have been 
developed 
RQ a\_4. Nora (8 On £ < OW 
Pla9s)=1-e (3) 1, Ce ae 
42 =0 
t me 
+5) VY } 
all 0. 2 Q 
ac = OO! > Ww 
=“. S > 3 Ty, (%2) foun 
de=1 


ememeetatiVely accUrace approximations to true hit proba- 
Pritles using Only the first five terms of the Series. 
Greater accuracy can be achieved by considering more terms. 
To have an accurate approximation to such a hit probability 
is beneficial to both weapon systems analysts and military 
commanders, because of the difficulties involved in comput- 
ing actual hit probabilities. The above equations yield a 
good approximation to hit probabilities using simple mathe- 
See ideameable Or modified Bessel functions, which is not 


titi reulle to obtain. 
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APPENDIX A 


COMPARISON OF HIT PROBABILITIES 


a ee 


ioe .26/ .267 w i ss i # # * # 
2.0 5 We Oe 40 39 7 a ey ‘2 a # # 
Pree 2508-2565 ./788 .786 .498 .433 ee a % 
Om eo eGo) 6.9660 °° 6.614 .803 .578 .450 sf si 


ee eee e998 4970 .969 .835 .813 .643 .460 


Pace eagles tor R Less Than a 


a iC ZO S20 4.0 Sy5 0) 
\3 app act ap Deeace appa act app ace App act 
i 


te ce er i et ep ee rer ene ee 


IenOne 2267 .267 = % % % % % % + 
wow .082 .082 .393 .397 a 2 % % # 4 
emOee O11 .011 .113 .113 .399 .433 : % % % 
PCemeOOC MOON. 2015 7.015 .121 .126 .366 .450 il * 


PCOOCm-OOOmm_OOL O01 2016 2.017 1.120 .133 .325 .460 


BDEguiem2—-s lable for R Greater Than a 


These tables compare the approximation to the actual hit 
probabilities (app) and the actual hit probabilities (act) 
themselves. Only the first five terms of the infinite series 
described on page 1] were used for the approximations. Accu- 
racy can be improved in all cases by considering more terms 


of the series. 


Us: 





BIBLIOGRAPHY 


Capos twa, Approximate Circular and Noncircular 
Gftsee Probabilities of Hitting," Operations Research, 
Valo pee 7 L-6255 1964, 


Eckler, A. R., "A Survey of Coverage Proglems Associated 
with Point and Area Targets," Technometrics, v. II, 
NOewo me Pee lDOl—580, August 1969. 


Guenther, W. C. and Terragno, P. J., "A Review of the 
Literature on a Class of Coverage Problems," Annals 
SecmdioniaricaleStatlstles, vs 35, pp. 232-260, 1964. 


Rand Corporation, "Table of Q Functions," Research Memo 
BaMauneNOn.n fJ57,008 Janvary 1950. 


CimiibbanGde bye ce, “Integral of the Bivariate Normal Dis- 
PeroMetonwoves an Offset Circle," Journal of the 
AiCm@mmcaneotatdistical ASSOCiation, Vv. 5/7, No. 300, 
December 1962. 


Luke, Y., Integrals of Bessel Functions, McGraw-Hill, 
OIG 2. 


ity) 





TP eAaheDisTRIBULICN LIST 


Defense Documentation Center 
Cameron Station 
Alexandria, Virginia 22314 


Library, Code 0212 
Naval Postgraduate School 
Menterey, California 93940 


Professor James G. Taylor, Code 55 Ta 
Department of Operations Analysis 
Naval Postgraduate School 

Monterey, California 93940 


Department of Operations Analysis (Code 55) 
Naval Postgraduate School 
Monterey, California 93940 


Major George E. Wrockloff III, USA (student) 


yoode blossom Drive 
Sam Antonio, Texas 78217 


20 


NO: 


Copies 





Unclassified 
Security Classification = 4, 
DOCUMENT CONTROL DATA-R&D 


(Security classification of title, body of abstract and indexing annotation must be entered when the overall report is classified) 





1 ORIGINATING ACTIVITY (Corporate author) 2@. REFORT SECURITY CLASSIFICATION 


Naval Postgraduate School Unclassified 


3 REPORT TITLE 


Pee ROXIMATION TO OFFSET CIRCLE PROBABILITIES 


4 OESCRIPTIVE NOTES (Type of report and inclusive dates) 


fieeeer's Thesis (March 1971) 


. AUTHOR(S) (First name, middle initial, last name) 


George Edmund Wrockloff, III 
Major, United States Army 














76. NO. OF REFS 


6 






8a. CONTRACT OR GRANT NO. 9@. ORIGINATOR'S REPORT NUMBER(S) 


b. PROJECT NO. 


c. 9b. OTKRER REPORT NO(S) (Any other numbers thal may be essigned 
this report) 


10. DISTRIBUTION STATEMENT 


Mepeoved for public release; distribution ia eaemaetee al ¢ 


»- SUPPLEMENTARY NOTES 12. SPONSORING MILITARY ACTIVITY 


Naval Postgraduate School 
MOMEe wey mecalirornta, 9 3940 





13. ABSTRACT 


If a weapon system is fired at a circular target and the impact 
Peiieeroedastitiouted as a bivariate normal random variable, it is not 
Pomeereubarly dirficult to determine hit probability if the expected 
impacc point is at target center. If, for some reason, the expected 
mmpact point is offset, the problem of determining hit probabilities 
becomes quite complex. Herein is developed a method to approximate 
ewen a hit probability. 


Se PL ee rT eh eT La fae JI zm Ser re ee eee 


DD werctore ntl Cade St ol 91 Unclassified 


S/N 0101-807-681} Security Classification canon 





Unie Vassi tied 


Security Classificetion 





LINK A —— 
A ois ee ee 


Oercet Circle Probability 
Pooroximation to Hit Probability 


‘Bessel Function Approximation 





— = i - 
DD (2e".1473 (sack) 
S/N 0101-+807-682}] 22 Security Classification A-31409 














Thesis 


126479 





W925 Wrockloff 

ee | An SPP roximation to 
Offset Circle Probabj- 
lities. 


¢ 


t 


©6970 





thesW925 


DUDLEY KNOX LIBRARY 





